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Abstract. 

We present exact results for the spectrum of the N*'' power of the Laplacian in a 
bounded domain. We begin with the one dimensional case and show that the whole 
spectrum can be obtained in the limit of large N. We also show that it is a useful 
numerical approach vahd for any N. Finally, we discuss implications of this work and 
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1. Introduction and main results 

Recently there has been a growing interest in the problem of finding the spectrum of 
large powers of the Laplacian in bounded domains [H [21 [3]. In one dimension, which 
is the case that will be presented in detail here, the problem is simply that of finding 
eigenfunctions and eigenvalues to the equation 



where u^''\x) is the k^^ derivative of u{x). Due to the invariance of (fT|2|) under 
the symmetry x —x, the eigenfunctions are characterized by a definite parity, i.e. 
u{x) = ±u{—x). In order to simplify the discussion it is useful to consider separately 
the eigensystem of ( fT|2l) for the even and odd eigenfunctions, which we will call from 
now on the even/odd spectrum. 




(1) 




(2) 



The spectrum of powers of the Laplacian in bounded domains 



2 



From a mathematical point of view, Ref. [Tj has shown that the determination of 
the spectrum of can be related to four different problems - the spectrum of certain 
positive definite Toeplitz Matrices, the norm of the Green kernels of A^, the best 
constants in certain Wirtinger-Sobolev inequalities, and the conditioning of a special 
least squares problem. From a physical point of view, the interest in the spectrum of A^ 
comes from many directions. Classical problems such as diffusion and wave propagation 
require knowledge of the spectrum of A (i.e., = 1), which is also related to recent 
problems such as diffusion limited aggregation and chaos. Problems in elasticity theory 
often deal with A^ [4j, and the fractional Laplacian (i.e., when N is not necessarily 
an integer) appears naturally in stochastic interfaces and Levy flights [HI E] and so 
knowledge of the spectrum allows progress in the understanding of anomalous diffusion 
and first-passage problems of a Levy flyer. 

The case = 1, emanating from classical diffusion problems, is trivial. The 
even spectrum is composed of the set of eigenvalues {[vr(j + 1/2)]^^ and their 
corresponding eigenfunctions cos(7r(j + l/2)x) (or (vrj)^ and sin(7rjx) respectively, for the 
odd spectrum). The case N = 2, which is encountered when analyzing the vibrational 
modes of a clamped rod, is also rather simple and discussed in [21 [3]. The even 
eigenfunctions are of the form u{x) = cos(Ax) + sin A cosh(Ax)/ sinhA, where A is any 
solution of the transcendental equation tan A + tanhA = 0. The odd eigenfunctions 
are of the form u{x) = sin(Aa;) — sin A sinh(Aa;)/ sinhA, where A is any solution of 
tan A — tanhA = 0. The allowed values of A cannot be written in closed form, 
although they are easily studied numerically and asymptotically. For = 3, a special 
symmetry simplifies the problem [21 [5]. As a result, the eigenvalues of the whole (odd 
and even) spectrum are exactly X2j = (jtt)^ and A2j+i = ((j + l/2)7r + 6j)^ with 
5j ~ 4(_i)i+ie"('^^/2)(2j+i) ig^^gg J |2]. Finally, for A^ > 3 it was shown that 

the eigenvalues of the whole spectrum are asymptotically known [2j. More precisely, 
there exists a jo such that for j > Jq 



A, = I ^ {j + 1) + - # cos [n {j + i) cos f ] + O (e^- #) | , (3) 

for A^ even, while for odd A^'s a slightly different expression exists 

{N 2N 
^J + ^e---tsin[vr,cos^]+o(e---lf)| . (4) 

However jo can be in practice rather large and thus the lower part of the spectrum, 
which is often the most interesting in physical problems, cannot be explored using ([3l 



^■+1 / . , 1 



2N 



Note also that Eqs. fl3|4j) do not provide any information about the eigenfunctions. 

Recently, it was shown [H [21 E] that for large A^ the eigenfunction with the 
lowest eigenvalue (so called "the ground state" borrowing the terminology of quantum 
mechanics) is given by Vq{x) = (1 — x^) and its corresponding eigenvalue is Ao ~ 
\/2(2N)\ with corrections which are of order 1/A^ compared to the leading order term. 
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Actually, in |3j many more terms were calculated for the lowest eigenvalue, namely 



3 25 375 8197 ^ / 1 



" ^ ^ 16N 512A^2 8192A^3 524288iV4 V A^^ 

Also, it was shown in [3] that the first excited eigenvalue is 



(5) 



(6) 



The present situation calls for more results concerning the spectrum of A^. In this 
paper, we show that in the large A^ limit, the eigenfunction of with absorbing BC 
can be written as certain associated Legendre polynomials. This also suggests, and is 
shown to be true below, that this basis is a good choice for diagonalizing the differential 
operator A^ also when A^ is not large. To be more precise, if we consider the normalized 
associated Legendre polynomials 



we show that the eigenfunctions of even spectrum are just given by 

Vjix) = P^^^,^ (x) + O (^-^^ for J = 0,1,2,.... (8) 

This result has a simple origin. In order to get some insight into it, recall that the 
associated Legendre polynomials form an orthogonal basis for functions in the interval 
[— 1, 1], that is 

P2N+^ix)P2N+Jix)dx = 6,, J = 0, 1, 2, . . . , (9) 

-1 

and they all obey the BC ([2]). This means that without loss of generality any 
eigenfunction of the aforementioned problem can be written as a linear combination 
of them. In particular, any normalized even function can be expanded using only 
even eigenfunction as follows (similarly, for the odd spectrum one needs the odd 
eigenfunction) 

oo 

U{X) = J2U^^^P2N+2AX), (10) 
3=0 
oo 2 

with ||f/(x)|| = Yl [U^^^ = 1. It is also well known [3 [3] that the ground state can 

j=0 

be obtained from a variational approach, namely it is obtained as the minimum of the 
functional 

/ U (x) [A^U (x)] dx 

X{U} ^ (11) 

J[U{x)fdx 
-1 

over all functions U{x) that obey the BC ([2]). In [H [3], it was proved that the ground 
state in the large A^ limit is simply given by Vo{x) = (1 — x^)'^, which is just p2^{x) 
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up to a numerical constant. This result is already consistent with Eq. ([S]) for j = 0. 
The next step is to look for the first excited state, which can also be obtained from the 
variational approach, as it is just the minimum of X{U} as given by Eq. (fTTi) over all 
functions U{x) that obey the BC ([2]) and are orthogonal to Vq{x). This means that in 
full generality, vi{x) can be expanded in the form 

oo oo 2 

V,{X) = J24'''P2N+2A^) ^ith Y.i'^?] =1. (12) 

i=i j=i 
in the large N limit, since the set < P2N+2j (^) f spans the (even) subspace which 

I ) j=l...oo 

is orthogonal to vq{x). 

Before providing a full proof below, we present a heuristic argument which can give 
some insight into the result ([8]). We claim that every zero crossing of a function U{x) 
that obeys the BC ^ adds a factor of order to the functional A{f/} (11 II) due to the 
presence of the large derivative, for large N . Recall that P2N+k exactly k roots in the 
interval [—1, 1] [8j. This means that in the general expansion ( |T2l) . the best candidate 
to provide a minimum for A{?7} is simply wj^^ = 1 and f^-^^ = 0(1/A^), as it has the 
minimal number of roots in the interval. Finally, this heuristic argument that yileds 
vi{x) can be recursively generalized to any fj(x), and the result (IHl) is obtained. 

In the rest of the paper, we first provide a proof of the statement ([8]). We use this 
result to obtain some explicit expressions for the eigenf unctions and eigenvalues in the 
large limit. Then, we will show that expressing in the basis of the associated 
Legendre polynomials does not only diagonalize it for N oo, but is also a very useful 
basis for numerically evaluating its spectrum for any A^. Finally, we discuss possible 
extensions and implications of this work. 

2. The Spectrum in the large N limit 

In order to prove the statement (IHl), we start by expressing the matrix elements of 
the operator A^ in the basis of the associated Legendre polynomials ([7]) 

1 

\x, (13) 



XN ^ / p2N / X 



-1 



aN p2N / X 
^ ^2N+2j l-'-J 



After some algebraic manipulations which are summarized in Appendix A, one finds 

XN ^ / (2i)!(4jV+2m)! V(4jV+4m+l)(4jV+4j+l)(-l)^+^- ^ (-4)'(2jV+2»)!r(2jV+j+»+|)r(»+|) 
m,j V (2m)!(47V+2i)! 2^^+^r(N+h) h \(2i)\f(i-i)\r( N+i+'^] 

C14) 



^-J,-N,2N + ^+J + i \ 2N + m + i-m,N+l 

' i + l,i + l ) \ 2N+l,N + i + l 

where T{u) is Euler's Gamma function, and 3F2 ( ; x ) is the generalized 

hypergeometric function [S]. Note that the matrix A^ is symmetric (A^^ = A^^), 
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as can be shown by integrating by parts fll3p N times and keeping in mind that the 
associated Legendre polynomials obey the BC ([2]). 

In the following, we will need two results regarding this infinite dimensional matrix. 
First, we show in Appendix B that the matrix elements behave asymptotically for large 
as 



A 



N 



'(2m)! 



2^- 



(2j)! (m-j)!(2j)! 
for m > j. Retaining one more term for the diagonal elements gives 



-lfV2{2N) 



(AN) 



2i 



^ 3 + 4j + 8j2 

1 — + O 

16N 



1 + 



1 

iV2 



V ^ ' V - (2j)! 
This means that to leading order the matrix has the following structure 

/ 1 J_ /i 



(15) 



(16) 



A 



TV 



V2{2N)\ 



1 

V2 



V(2j)! 

2ij\ 23-''/^{j-iy. 



32jY4 
3 



2j-ii/233/2Q-„2)! 



V(2i)! 



2j-7/2Q_l)!^^ 



2j-ii/233/2(j_2)!^^ 



2.1. Proof for the diagonalization 

We will now prove by induction that in the large A^ limit the spectrum of this matrix 
is just its diagonal elements so that it is essentially diagonal in this basis to order 1/A^ 
and thanks to some extra structure, this is true actually up to order l/N"^. 

Let us suppose that this statement is true for the submatrix of size K defined as 
A^f = \ • \ and then we will show that it holds for A^'^"*"^ as well. Take 

I J m,j=0,--,K 

the trace of the finite [K + 1) x [K + 1) submatrix. Using Eq. flTB]) it is given explicitly 
by 



^ (2j)! ^ ^ V^(2ivj. ^2;^ ^2)! 



TrA^'^+i = (-1)^ v^(2Ar)! ^-^^ = {-if V2{2N) 

j=0 

This simply means that Af^'^{^^^^ — TrA^'^+-^ = O {^)- Note also that A^_,_;^^^_,_]^ is 
the only term in the submatrix A^'^+^ that is of order (2A^)!A^^'^"'"^. This together with 
the fact that the trace is just the sum of all the eigenvalues, means that A^^^ is an 
eigenvalue of A^'^+^, and it is the maximal one. It is easy to check that the column 
vector vk+1 = (0, ■ ■ ■ , 0, 1)"^ is the eigenvector that corresponds to this eigenvalue, up 
to corrections of order O {l/N"^), since 

/\N,K+^ _l n^,,/^ro,Ani ( v/(2X+2)(2X+l)(4jV)^^ (4^)2(K+i) \^ 

^ Vk+1 - l-J-J (,^iVj! I 2^+i{i<r+l)!' ■ ■ ■' 2(2K)\ ' (2E'+2)! ) 
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(19) 



1 

We have succeeded in finding one eigenvalue and one eigenvector, and in order to 
complete the diagonalization we need to consider the upper left K x K submatrix. 
However, this is nothing but A^''^ which is by assumptions diagonal in this basis. 
QED. 

2.2. Higher order results 

Note that all the arguments of the proof above are valid up to corrections of order 
O {l/N"^). This implies that the diagonal elements of are indeed the eigenvalues up 
to O {l/N"^). This fact becomes useful when considering higher order corrections, since 
we know that Eq. (fT6ll already gives the first subleading corrections, without any need 
to diagonalize the matrix. 

The same reasoning and a quick look at Eq. f|T9|) lead to the conclusion that in 
order to get the first correction for the eigenfunction (which is of order l/N"^) we need to 
diagonalize blocks of 3 x 3 (with the exception of the ground state where the 2x2 upper- 
left submatrix suffices). Thanks to the special structure of the asymptotic behavior of 
the matrix, it is sufficient to know the ratio K^^^j^^/ K^j^^^^j^-^ up to order 1/A^, namely 



A«;,,« 32 V »N ) \N-' 



and so we get 

Vj (x) = P2N+2j (^) 



+ ^ { m (2j - l)]'/'P|^+2,_2 (x) - [{2j + 2) (2j + l)]'/^P|^+2,+2 {x)}+0^^ 

In turn this allows to know the next two subleading orders for the eigenvalues, that 
is up to order After some algebra we get 



Xj = V2{4:N) 



2i 



3 + 4j + 8f 75 + 3Uj + 672 f + 832 f + 192 f 



.{mi 

{2j)\ 16N ' 1536iV2 



+ 



1125 + 5028j + 11960j2 + 5120j3 + U16f - 5888j5 - 5l2f / 1 
^ 24576iV3 ^ V iV4 



.(22) 



Repeating the same method, we can also obtain even higher order terms for specific 
states without serious difficulties. For example, the first two eigenvalues are given by 



Ao = V2{2N) 
Ai = 8V2N^ {2N) 



^ 3 25 375 8197 694941 / 1 

" ION ^ 512A^2 + 8192A^3 ^ 524288iV4 ~ 8388608iV5 ^ [m 



(2 



^ 15 705 7083 _ 963989 _ 68014113 ^ / 1 
~ Ton ^ 512A^2 + 8i92Ar3 ~ 524288iV4 ~ 8388608iV5 ^ 

which are consistent with the results of [3], where these quantities were calculated up 
to for the ground state, and up to 1/A^ for the first excited state. 

Recall that these results are exact when — oo. However, for a finite A^, there 
always exists a certain j for which these first terms are not sufficient, and in order to 
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obtain the eigenvalues/eigenvectors one needs to know more terms of the expansion. 
This should have been expected since as has been conjectured in [3] and proved in [2] 
there always exists a jo (which can become very large actually) such that for j > jo the 
eigenvalues approach asymptotically (| + vrj) for odd and + vrj) for even 
N. Such a transition to "simple" eigenvalues necessarily implies that many terms in the 
large N series should be retained as we get to higher excited states. 

So far, the results mentioned above apply only to the even sector of the spectrum. It 
turns out that the odd part can be obtained by simply transforming j ^ j + | in all the 
formulas for the even part, namely Eqs. ([8]) and (fT5|) - fl22|) . As a consequence the proof 
for the even part can be simply repeated for the odd part. Therefore, we essentially 
fully cover the two sectors of the spectrum for large A^. An interesting conclusion, that 
was conjectured already in the introduction and is now proven, is that in the spectrum 
each extra zero-crossing of the eigenfunction "costs" a factor A^ in its corresponding 
eigenvalue in the large- A^ limit. 



3. Numerical approach 

The proof presented in the previous section shows that the orthonormal set ^p2N+2j 
diagonalizes the operator for large A^. However, from a more general point of view, 
this orthonormal set is a complete basis for even functions that obey the BC ([2]) and 
therefore, by representing A^ in this basis we loose no information. For that reason, a 
direct diagonalization of the matrix A^ ^ should provide the eigensystem of A^ for any 
A^, even though it is not clear a priori if this approach is efficient. On the other hand, 
since we proved that this basis diagonalizes A^ in the large A^ limit, we expect it to be 
a good starting point for diagonalization for any A^. Indeed, we now show that one can 
get very accurate results for finite A^ by numerically treating the matrix A^^- f|T3|) . 

In order to show the high rate of convergence we can consider the case A^ = 1 
directly. Note that for A^ = 1 we know exactly the whole spectrum, which is 
|cos (l + nvr) ^ with its corresponding eigenvalues + nvr) . The following 
expansion holds 

oo 

COS (f + nn)=J2 CiPi+2^{x) , (25) 

with 



i=0 



(2^ + 4)! 



(4^ + 5){2(2n + l)7rj^, 



(2n+l)7r 



i+2 I 2 



+ 



+ 



2 + -M- + \f] 32.+2 (^) + vr^ (n + i)^-,2 (^) } 



(26) 



' 22^+2 (22 + 1)'^+=^ 

The coefficients Cj decay faster than exponentially to zero for large which means that 
one practically needs a small number of coefficients in order to get a good approximation 
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of the real eigenfunctions. Note that when diagonahzing the matrix Aj ^ ((Eq- fll3p with 
= 1), we exactly calculate the coefficients Cj, and therefore it is enough to consider 
just a finite submatrix approximation to the full matrix Aj^m in order to get a good result 
for the eigensystem. A simple thumb rule of taking a finite submatrix, diagonahzing 
it and then keeping only the ffist half eigenvalues yields result with good accuracy. In 
Fig. [1] we present two examples where we evaluate the accuracy of this approach for the 
eigenvalues for = 1 and A^ = 50. In both cases we assess the accuracy by comparing 
the eigenvalues obtained from a truncated 51x51 matrix, i.e. A^^^, to those obtained 
from A^^°. For A^ = 1, since we know the exact eigenvalues, namely \j = [n (j + |)] ^, 
we can also have an absolute check on the accuracy. As one can see the simple thumb 
rule of keeping only the first half eigenvalues amounts to a relative error of at least 
~ 10~^°. Keeping in mind that the eigenvalues grow very fast with A^ and j (here we 
only look at ratios in order to have a better evaluation) this is reassuring. 



~. — ; 1 1 r- 



-(a) 



CO 
X 

o 
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10 
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Figure 1. A plot of the estimated accuracy for the first 50 eigenvalues obtained from 
a truncated 51 x 51 matrix, i.e. A^-^^, for (a) = 1 and (b) N — 50. In both cases 
we assess the accuracy by comparing to the first 50 eigenvalues of A^^°. For = 1 

we also compare to the exact result Xj = [tt (j + 5)]^- Note that the two estimations 
agree quite well with each other. 



Concerning the eigenfunctions, the situation is of course more delicate. In order 
to get sufficient convergence throughout the interval [—1,1], we had to keep around 
1/3 of the states. In Fig. [21 we compare the lO*'^ excited state Vio{x) with the exact 
result cos(^x) for A^ = 1. As can be seen, the two functions agree very well and the 
difference reveals a discrepancy of order 10~^, which is mainly concentrated close to the 
boundaries. 

The latter discussion has an important implication. If the above procedure 
converges rapidly for A^ = 1 then it certainly does so for higher values of A^, which 
means that this numerical procedure can be useful in general. Having established the 
fact that the numerical scheme yields results with controllable accuracy, we can use the 
outcome in order to assess the various analytical results obtained so far. In Table [H we 
compare the numerical results for the eigenvalues obtained from the truncated matrix 
of 50 X 50 (i.e., A^^°) to the large A^ 4*'^ order formula we obtained in Eq. (l22l) and 
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Figure 2. A comparison between the 10*'' (even) excited state obtained from the 
numerical scheme using a truncated matrix of 50 x 50, i.e. A^^'^^, denoted here viq{x) 
and the exact eigenfunction cos{^^x) for the case = 1. In (a) we plot the two 
functions, where they seem to agree very well, while in (b) we plot the difference of 
the two in order to allow for a more quantitative evaluation. Most of the error is 
concentrated near the boundaries, where it is of the order of 10~^. The integral of the 
absolute value of the difference is J^^ \vio{x) — cos{^^x) \ dx — 3.43 x 10^^. 



j 


A5'""V%/2(2A)! 


^num 1 yAth 


A™"V [7r(i + 3/4)]^^ 


1 


0.99627 


1. 


7.90282 X 1012" 


3 


6.33108 X 10^ 


0.999996 


1.87819 X 10^2 


5 


5.35907 X 10" 


1.00007 


2.91917 X lO^-i 


7 


6.04205 X 10^8 


1.00309 


1.80442 X 10'*'* 


9 


1.7643 X 10^3 


1.0346 


2.82683 X 10^7 


11 


1.86415 X 10^7 


1.27419 


3.42899 X 10^2 


13 


8.751 X 10^0 


-42.7013 


6.25995 X 1028 


15 


2.10073 X 10^4 


-0.337001 


7.05830 X 1025 


20 


5.86444 X 10^1 


-0.0158765 


4.14045 X 1020 


25 


1.31862 X 10^8 


-0.00182839 


1.47275 X 10*^ 



Table 1. A comparison between the different results for the eigenvalues of the even 
spectrum for A^ — 50. A""™ that results from the numerical procedure by using a 
truncated matrix of 50 x 50, A*"* that come from Eq. ([22|) and the asymptotic result 
[7r(j + 3/4)]^^ [2]. As can be seen, the first 11-12 eigenvalues are reasonably captured 
by A**'', while the asymptotic formula is still far from giving a good estimate even for 
j = 25. The negative values reflect the fact that A**'* is no longer a good approximation 
for j > 12. 



to the asymptotic result [7r(j + 3/4)] of [2J for N = 50. First, notice how fast the 
eigenvalues grow already for = 50. Second, one can see the large N formula (!22l) 
describes well the low states and starts deviating around j ~ 12 where they become 
negative. Table [1] also shows that when A^ becomes large, the critical jo above which 
the asymptotic result of [2| is valid becomes large itself. 
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4. Summary and discussion 

In this paper we calculate the spectrum of the N^^ power of the Laplace operator 
in the large limit. We obtain explicit expressions for all the eigenfunctions up to 
order and for all the eigenvalues up to order for both the even and the 

odd parts of the spectrum. In addition we show the general structure of the series 
which provides corrections to the eigenvectors of order l/N"^^ by diagonalizing blocks 
of (2m + 1) X (2m + 1) for large A^, and keeping consistently all term up order 
in the matrix elements that contribute. The expansion we obtain is consistent with 
some previous known results for the ground state [H |3] and for the first excited state 
[3]. In addition, we show that this approach can be useful for small A^'s, where a direct 
diagonalization of the matrix representation of the operator in the basis of the 
associated Legendre polynomials converges very fast to the exactly known spectrum for 
N = 1. This of course implies even faster convergence rate for higher values of N . 

Some applications of this result comes to mind. For example, in the context of 
the Wirtinger-Sobolev inequalities, which are the generalizations of the string {N = 1) 
and rod (A^ = 2) inequalities [9], one can benefit from the results derived here. A 
Wirtinger-Sobolev inequality is simply an inequality of the form 
11 1 

J [U{x)fdx <C J U{x) [{-AfU{x)] dx = C J [U^^\x)Ydx, (27) 
-1 -1 -1 

where U{x) is required to satisfy certain additional conditions, such as the BC above ([2]). 
The variational formulation ( fTTi) shows that the smallest eigenvalue Aq of the operator 

indeed gives the best constant C, with C = 1/Xq. However, knowing the whole 
spectrum of A^ gives much more information. For example, suppose that we know 
that a given function U{x) is even (odd) and that it crosses zero 2j (2j + 1) times. In 
that case, the best constants in the Wirtinger-Sobolev inequalities can be improved and 
become C = (or C = l/Xj^i for the odd case). For large A^, these quantities are 
determined using Eq. (!22|) . and for any A^ they can be calculated efficiently using the 
numerical procedure described above. The key result here is that each zero crossing of 
U{x) in the interval lowers the constant in the inequality fl271) by a factor proportional 
to A^ (when A^ is large). 

Interesting generalizations of the present results include fractional values of A^ OE], 
other boundary conditions and more than one-dimensional Laplacian [3J. Already for 
large, but not necessarily integer A^, the results we derived for the eigenvalues hold even 
for the fractional case. Moreover, we believe that the approach presented here can be 
adapted to yield the whole spectrum of those operators in the large A^ limit, and through 
a rapidly convergent numerical approach for any A^. The main difficulty is to find a set 
of orthornormal polynomials for the case of fractional Laplacian which plays the same 
role as the associated Legendre polynomials for integer A^'s. 

A more challenging class of problems are those where the domain is not simple, such 
as a rectangle or a billiard. An approach in the spirit of the present one might allow 
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some analytical insight through a expansion into problems that are traditionally 
heavily numerical. For example, for the 3D Laplacian case the problem can be posed 
as solving [3] 

- ^^^^^)\ru{r)) = i-lfXiruir)) , (28) 

with u{r) satisfying the BC ([2]). One can show that the appropriate set of orthonormal 
polynomials are 



v,{r) oc r^(l - ^2)iVp{27V,L+i/2)^2r^ - 1) , (29) 



where P^°''^\x) are the Jacobi polynomials. Using this observation, one can repeat the 



same procedure as presented here in order to obtain the whole set of eigenvalues. 
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Appendix A. Derivation of the matrix elements A 



N 



In this appendix we present how to obtain the matrix elements A^^- given in Eq. f[T3]) . 
The starting point is the explicit series representation of the associated Legendre 
polynomials that can be obtained for example from Equation 8.812 in Ref. [8j 

^2iv+2, W - 22JV+2, X) Z^(2£)!(j-£)!(2iV + j+£)! ' ^ 
We can now explicitly act on this representation with the operator 

^ nN+2j[X)- ^2N-,2j 2-. ^2£)\ij-i)\{2N+j+£)\ ^ ' 



'-iy~^ (4iV + 2j +2£)! 
!£)!(j-£)!(2iV + j+£) 

(A AT ^ OA , new ^^^W) .N~m 



1 ^ i-iy-' {AN + 2j + 2i)\ l^ 'N\{-iy'"'^{2N + 2i-2my. ^ 



22^+2i (2£) ! (j - tj ! {2N + j + tj\ m\ {N - my. (2£ - 2m) ! 

This in turn can be written as A^P2^jv+2i (^) ~ Tl, with 

1=0 

(,,^) ^ (-l)^(2iV + 20! ^ {-iy^\AN + 2j + 2m + 2iy A^! 

i 22^+2i(2i)! ^^{2m + 2iy.{i -m-iy.{2N + i +m + iym\{N -my.'^ '' 



Using the integral 7.126-2 from Ref. |8j 

^ ' 24'vr(7v+i)r(7v+i+|) (2^-2^)!'^ 2N + l,N + i + l ' ''^ ' 

we can simply obtain Eq. (IT^ . 
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Appendix B. Large N asymptotic expansion of the matrix elements 



In this appendix we show how to obtain the large expansion of the matrix elements 
A^^- summarized in Eq. f|T5l) . In order to facilitate the calculation we assume from now 
on that m > j, since the matrix is symmetric and so the case m < j is trivially obtained. 
The first step is to simplify the hypergeometric function 3F2. We use formulas 7.527-1 
and 3.259-3 from Ref. [8] to obtain the following integral representations 

2N + j + l-j,N + l 



3-^2 



2N + l,N + i + l 



;1 



7.527-1 



-2 e-(^+i)-2Fi 



2j 

1 _ 1 



2N + J + 1, 
2N+1 



dx 



(B.l) 



J {l-e—y~' 2 e-(iV+l)- / 2 {l+y)-(^^+^+'') [l+{l-e-)yy dydx 



3.259-3 



B(jV+l,i+l)B(2Af+j + i i-j 



where B{x,y) is the Beta function. In addition we expand [1 + (1 — e ^)yy 



E 

e=o 



'A(j-ey. 



[l — e ^) and so 



3-^2 



2N + j + l-j,N+l 
2N +1,N + 



E 



;1 



1 

)'+^-2e-(iV+l)^dx 



FT" 



After using J y°- [1 + y)^ dy = 


get 

" 2A + j + i,-j,A+l 
2A+l,A + z + f 



B{N+lMi2)B[2N+j+^,^^j'l 

00 

r(a+i)r(-i-a-6) J _ g_,^a ^.fo,^^ = 5 (a + 1, 6) we 



;1 



1 1_ . 
'2'2 



B(Af+l,i+F) 



(B.3) 



Last, for large N we can expand the Beta functions and get to leading order (subleading 
orders are not difficult to obtain) 



(B.2) 



3-^2 



2A + j + i,-j,A+l 
2A+l,A + i + | 



\/7r2-^i 



;1 



* + 2' -3 

1/2 



;2 



2iAJT(i-j) 



O 



Next, we expand the factorials in ft-"''^'' as defined in Eq. (1A.3P 



TT- 



(2i)!j! (j - 0!r (j + i) 



SO that 



(4jV+4Tn+l)(4jV+4j+l)(4iV+2m)!(2j)! 



r i 



(4Af+2i)!{2m) 
72 (2A)! 



24iV+ir('7V+i) r(Ar+i+| 



-^^(_l)JV+J+«22» + 2m^2i + mr^j+l^ 

(2*)!j!0-i)!r(i+i) 



(B.4) 



(B.5) 



(B.6) 



After combining all these results, and some simple additional algebra we obtain Eq. fll5p . 
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